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Abstract 

A lattice QCD study of the strong decay width and coupling constant of decuplet baryons to an 
octet baryon - pion state is presented. The transfer matrix method is used to obtain the overlap 
of lattice states with decuplet baryon quantum numbers on the one hand and octet baryon-pion 
quantum numbers on the other as an approximation to the matrix element of the corresponding 
transition. By making use of leading order effective field theory, the coupling constants, as well as 
the widths for the various decay channels are determined. The transitions studied are A —>■ n N, 
£* —> A 7r, Y,* —> £ 7T and £* —> Ett. We obtain results for two ensembles of Nf = 2 + 1 dynamical 
fermion configurations, one using domain wall valence quarks on a staggered sea at a pion mass of 
350 MeV and a box size of 3.4 frn and a second one using domain wall sea and valence quarks at 
pion mass 180 MeV and box size 4.5 fm. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70.Dj 
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I. INTRODUCTION 


The study of resonances from first principles using lattice Quantum Chromodynamics 
(QCD), has progressed significantly. Most of these studies are based on the Liischer approach 
DEI and extension thereof [SHU that extract scattering lengths and phase shifts from 
discrete energy levels in a finite volume. The approach has been generalized to the case 
of coupled channels [SHE] and three identical boson scattering [18], and a growing number 
of studies is being carried out in the meson sector. A pioneering study of meson-baryon 
and baryon-baryon scattering lengths was already carried out twenty years ago [19] and 
more recent studies include those by members of the NPLQCD |2_0. 2T] and HALQCD [22] 
collaborations and other groups [22]. Despite this progress, the application of the Liischer 
approach to baryon resonances has been limited since the method requires very precise data 
for multiple spatial volumes or various reference frames of different total linear momentum 
making it computationally very demanding. 

Another method to study hadronic resonant decays from lattice QCD was proposed in 
Refs. (2T[ '25] and successfully applied in the study of meson decays [26i Eg. A first applica¬ 
tion of this transfer matrix method to baryons was carried out in Refs. [28l [29] . The transfer 
matrix method as applied here allows to extract the width of a resonant hadronic decay, 
if the resonance width is small as compared to the resonant energy and well-isolated from 
other decay channels. In such a situation the method allows us to extract the width from 
one kinematic point and it thus provides currently a computationally feasible calculation of 
the width in the baryon sector. This calculation can be seen as a first attempt to compute 
the width of an unstable baryon that allows us to learn about two-particle interpolating 
fields in the baryon sector and the associated technicalities and gauge noise. 

In this approach, one considers a purely hadronic decay of a baryon B* to a two-particle 
state. In the cases considered in this work, the two-particle state will be a meson M and a 
baryon B , as illustrated diagrammatically in Fig. [lj. 

We associate a vertex with the tree-level transition graph in Fig. [l] and the strength of 
the interaction at the vertex is measured in terms of an effective coupling constant g^s- 
We define this coupling to coincide with the coupling that appears in the leading order 
continuum effective field theory for the interaction term of the hadronic fields M, B and 
B* in the effective Lagrangian. This will be made more explicit later on in connection with 
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FIG. 1. Diagram for the transition B* «->■ MB 


Eq. 


( 10 ). 


In order to study a decay B* —>■ M B in the Euclidean quantum field theory we need 
to formulate it in terms of energies or hadronic matrix elements. In lattice QCD we use 
interpolating fields to create states with the quantum of the decuplet baryon B* and the 
octet baryon B and the meson M. We restrict our consideration to the two lowest-lying 
states with the desired quantum numbers, which we label by | B*) and | MB). If B* does 
not decay then its overlap with | MB) is zero and it is an asymptotic state of the theory. 
These states can then be thought off as the eigenstates of a non-interacting lattice transfer 
matrix T° defined by a Hamiltonian H 0 = \B*)(B*\ + | MB)[MB |. Our approach here is to 
study the overlap of the states created by the interpolating fields B* and M B for the case 
where the energy levels of these states are near-degenerate. The interaction Hamiltonian to 
leading order in the perturbation will then be given by H = Ho + \B*)(MB\ + | MB)(B*\. 

The off-diagonal elements of the Hamiltonian will be the overlap (. B * | H | MB). Thus 
our assumption is that a state | B*) created initially at time ti propagates in Euclidean time 
on the lattice to final time tf, makes one transition to | MB) at any intermediate time step 
t —> t + a on the lattice. If the (real valued) lattice transition amplitude is small in terms of 
the inverse propagation time (tf — ti) 1 and if the energy gap between the states | B*) and 
| MB) is sufficiently small then one can evaluate the overlap and relate it to the coupling 
constant and then to the decay width [25] . We stress that the propagator of the state B* 
is fully dressed and so is the propagator of the M B state including interactions between 
the two particles. A tree-level effective interaction Lagrangian can be written in terms of 
the fields B *, B and M [30] with coupling constant g^Bi which is related to the overlap 
(B* | H | MB) as will be discussed in section II. 


In order compute the overlap of these states we need to choose ensembles for which the 
energy gap 5 = Eb* — Emb is small in units of the inverse propagation time from initial to 
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final state 6 -C l/(tf — ti). Since on a finite lattice the allowed momenta are discretized the 
energies will not in general match. Thus, this condition will be only approximately satisfied 
for the ensembles we have at our disposal. A second condition that is required is that the 
propagation time tf — ti is sufficiently large compared to the energy difference between the 
ground state energy of B* and its first excited state as well as between the lowest energy of 
the B M system and its excited state with the same quantum numbers so that only the two 
lowest-lying states of interest dominate in the transition matrix element. The extraction 
of the overlap from lattice measurements is detailed in section II. The transition matrix 
element Ml(B* —y MB) oc ( B * \ H\MB) for the situation in which the energy levels of 
the two states are degenerate. Using Fermi’s Golden rule one can relate this decay matrix 
element to the decay width 


T^ b = 2-k\M{B'^MB)\ 2 P , 


( 1 ) 


where p is the density of states at the transition energy. As already mentioned, in this study 
we work to lowest order considering only a single transition amplitude and allowing for large 
enough time separation t/ — U so only the lowest states in the initial and final states give the 
dominating contribution. To this order we also neglect further elastic rescattering of MB 
in the final MB state. 

In our first study [28], we successfully applied this approach to study the A resonance 
using a hybrid action with domain wall valence quarks on a staggered sea. Here we extend 
our study to include the decuplet baryons E* and S*. In addition, we investigate the 
applicability of All Mode Averaging (AMA) j31j in improving the statistical accuracy using 
the A resonance as test case. In this work we also analyze an ensemble of domain wall 
fermions (DWF) corresponding to a pion mass of 180 MeVpJZj for which the energy matching, 
in particular for the A, is very well satisfied. The results based on this ensemble of Nf = 2 + 1 
DWF for the widths of the A, the E* and the S* as well as the results using the hybrid 
ensemble for the E* and S* constitute the first determination of the decay widths of these 
resonances using lattice QCD. 

The paper is organized as follows: In section [TT] we present the method and give the tech¬ 
nical details. In section m we show our lattice QCD results, in section [TV] we discuss these 
results and their relation to the decay widths and in section [V] we present our conclusions. 
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II. TECHNICAL DETAILS OF THE METHOD 


The method that we consider in this work was first described in Refs. [23, 25j [33J where 
it was applied to the study of meson decays. The method was extended for the case of 
the A resonance and first results were obtained using an ensemble of domain wall valence 
quarks on an Nf = 2 + 1 staggered sea, which we will refer to as hybrid approach [28] . 
In this section, we explain the technical steps involved paying particular attention to the 
description of the decays of decuplet baryons, which is the focus of this work. 


A. Lattice correlation function and normalization 

We consider the following strong decays of a decuplet baryon to a nreson-baryon final 
state: 


A —> 71N 

£* —» 7T £ , 7T A 

S*-MrS, (2) 


generically denoted by B* —> M B. Due to the isospin symmetry of the lattice action, we 
can choose any isospin channel for each case. In this study, we consider the A ++ , the E* + 
and the S* - with interpolating fields given by 




iPx 


j£ a ++(x) = e abc (u aT (x)C^u b (x)) u ac (x) 

J^ + (x) = A=e abc (u aT (x)C Yu b {x)) s ac (x) + 2 (s aT (x) CY u h {;x)) u ac (x) 

v 3 L 

J^-{x) = e abc (s aT (x)C Y d b (x)) s ac (x ), 


(3) 


where lower case Latin (Greek) letters denote color (spin) indices and C = i'y o 72 is the 
charge conjugation matrix. 

As interpolating fields for the meson-baryon states we take the product of the interpo¬ 
lating fields of the corresponding meson and baryon: 
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JM B (t,PM,PB) = Y, Mt, V) J%(t, 


y,z 


Jn+(y) = d(y)'y 5 u(y) 

Jn(z) = tabc (u aT (z) C 75 d b (z)) u ac (z ) 


i Z ) — ^/g € abc 


2 (u aT (z)C'y 5 d b (z)) s ac (z) 


+ (u aT (z)C^ 5 s b (z)) d ac (z) - (d aT (z) C 75 s b (z)) u ac (z) 


{u aT {z)C^s b (z)) d ac {z) + (d aT {z)Cj 5 s b (z)) u ac (z ) 


^s 0 (a) ^2 e “ bc 

J H o = e abc ( s Ta (z)C'f 5 u b (z )) s“ c (^). 


(4) 


Previous studies have shown that the two-hadron interpolating fields given in Eq. (JdJ) will 
primarily overlap with two-hadron states [31], while the single-hadron interpolating fields in 
Eq. (J3]) will have a dominant overlap with single-hadron states. 

In the following, we consider kinematics where the total momentum is zero, so in the Erst 
line of Eq. (J 3 J) we set P = 0 and Pm = ~Pb hi Eq. Q. Using these interpolating fields we 
build the two-point correlation functions for C b *-b *, Cb*-mb and Cmb-mb as follows: 


C B *-B*(tf-ti) — Tr |^-(l + 7o )Pk! 2 ^2( jl B*(tf,x) J B*(t i ,z))j (5) 

- U,q) = Tr ( 1(1 + 70 ) pT v (J‘ B .(t f ,x) e-W - 25 

( 6 ) 


x,y 


C M B-MB(tf - ti,kf,ki) = Tr I -(1 + 70 ) ^ {J B (t f ,x) J M {tf,y) Juit^y') J B {t h z)) 

\ S,y,y' 

x e ik f (y-x)-iki(y'-z) 


(7) 


with a fixed source location (ij, z). All correlators are defined to include a parity projection 
4 (1 ± 70 ). In addition, Cb*-b * and C b *-mb include a projector to spin 3/2, which at zero 
total momentum is given by 

Pfk 2 = $ik 1 - ^ 7*7 k , i,k = 1,2,3. ( 8 ) 
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To cancel the unknown overlaps of the interpolating fields with the states we construct 
the ratio 


R 


B* 

MB, k 


t f - t t , q, k f , hi) = 




Cb*~b*( tf —ti) x Cmb—mb( tf —ti,kf,ki ) 


(9) 


In this work we always consider the case kf = ki and |<f| = \kf\. The 2-point function 
Cmb-mb then only depends on \kf\ and the ratio can be characterized by a single vector 
k = q as RM B k = Rmb, k(tf — U, k) 


a. Alignment and polarization for B* — MB, momentum averages and angular momen¬ 
tum: Within leading order in effective field theory, a non-vanishing signal in q ) 

only arises when the relative momentum vector is aligned or anti-aligned with the spin pro¬ 
jection appearing in the correlation function, i.e. when q ■ e*. ^ 0, where edenotes the unit 
vector in the k— direction. The vertex for the fields B*, M, B in the effective Lagrangian 
following our notation is given by 

C,~gS B S;S,M“T"B (10) 

with matrices T“, which contain the Clebsch-Gordan coefficients for coupling isospin chan¬ 
nels. 


( T a )ik — {Ib* — 3/2, — i | Ib — 1/2, 1% — k; Im — 1 , I'm ~ a ) < 

i e {-3/2, -1/2, +1/2, +3/2}, k e {-1/2, +1/2}, a G {-1,0,1}. 

We perform our calculations with one unit of relative momentum, \q\ = 2ir/L, such 
that q = (±1,0,0)27 t/L or a permutation thereof and thus look at the six combinations 
C-b*~mb(^i ±27t/L efc), k — 1,2, 3. 


The B* — MB correlator is projected to its spin-3/2 component with P 3//2 . Moreover, 
the average over positive and negative momentum effectively means that for the MB state 
we use the interpolating field in its center-of-mass frame 


J 1 mb (t, P = 0, k = 27 r/ l)j = - ^2 P?/' 2 H mb (t, P = 0, ke^j - J mb (t, P = 0, -ke v 

i' 1.2.3 


( 11 ) 
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In a partial wave expansion, we find that the dominant state excited by the interpolating 
held given in Eq. ( [IT] ) will have orbital angular 1 = 1. The coupling of the l = 1 to the 
nucleon state with spin 1/2 are projected to the component with total angular momentum 
3/2. Thus the operators in Eq. (11) and the projected operators in Eq. ^ transform 
under the spin-3/2 representation of the Lorentz group in the continuum, which is subduced 
into irreducible representation A p = H + (positive parity) of the double cover 0 D of the 
octahedral group on the lattice (Table III in Ref. [35]). The irreducible representation H 
contains an overlap with higher partial waves. What simplifies the calculation at hand, is 
that we only consider the ground state at large Euclidean time, in which all channels but 
the desired one with l = 1, J p = 3/2 + are exponentially suppressed. 

All the standard spin-3/2 interpolating fields Jb* involve the spin structure (q T C 'fk q) q 
and the parity operation in the center-of-mass frame acts as 


P: 4.(t,P = 0)-> 7 o J|.(t,P = 0). 


( 12 ) 


The spin-3/2 projector, P 3 / 2 , and the projector to the component with definite parity com¬ 
mute, such that we have the trivial action 

P: i(l+ 7 o)W4'*(«.^ = 0)^i(l+7o)W4'*(«.^ = 0). (13) 

On the other hand, for the MB —state with the pseudoscalar meson held we have in the 
center-of-mass frame and with relative momentum k 

P : |(1 + 70) Jmb^, P = 0, k) ->• -i(l + 70) JmbP, P = 0, -k). ( 14 ) 

Since parity is a symmetry of the lattice action and the B* — B* and MB — MB two-point 
functions are even under parity, we expect the following relation to hold for the ratio 

R'MBp't ~~ til k) = R\lB(jt — til —k) ■ ( 15 ) 


b. Quark-connected and disconnected diagrams: The Wick contractions for the corre¬ 
lation function Cmb-mb can be represented by two types of diagrams, as shown in Fig. ( 2 ] : 
quark-disconnected (. D , upper right) and quark-connected (lower diagrams 61 , 2 , 3 ). 

We make two simplifications: i) We neglect the quark-connected diagrams Cj , C 2 and C 3 
and ii) approximate the quark-disconnected diagram by the product of expectation values 
of the individual meson and baryon propagators. This results in a significant reduction in 




FIG. 2. Diagrams representing different types of Wick contractions for Cb*-mb (diagram T) and 
Cmb-MB (upper right: quark-disconnected D , lower: quark-connected C\ ,2,3)- 

the computational cost. The contractions for diagrams T and D, as performed in this work, 
require one quark propagator and one sequential propagator through the source timeslice 
ti. In contrast, the full contractions for the quark-connected diagrams require all-to-all 
propagators. We thus set 


CMB-MB(tf — ti, kf, ki) —* CMB-MB^tf — ti, k]tf, k B ) 

= (CM-M(tf — ti, k M )) x Tr ^-(1 + 70 ) ( C'B~B(tf — ti, k B ))^j , 

( 16 ) 

where the meson and baryon propagators are Fourier transformed with independent mo¬ 
menta k^i and k B , respectively. In this form the M — B two-point function depends only 
on the squared momenta k 2 M and k 2 B and for Cmb-mb we only use combinations with equal 
modulus such that we can simply replace the dependence on the pair (fc M , k B )) by a single 
vector k. 

Thus, in the asymptotic region of large time separation (tf — t^/a 1, we can express 
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this approximation as 


— ti, kf, ki ) 


(f <x >> Z MB (k f , ki) e u) UMB (k f ) u M B(ki) + ... 

« Z M (k 2 ) Z S (P) e ~ E B(. t f~ t i) UB (k f ) u B {ki) + ... 


(17) 


B. Extraction of coupling and the matrix element 


For all baryons B* to M B we restrict the lattice Hilbert space to two states | B*) and 
| MB) as the dominant baryon and baryon-meson states, respectively [2Hj. In terms of this 
two-dimensional subspace, the transfer matrix is parametrized as 


t = e ~aE 



(18) 


In accordance with the assumption of a small energy gap we take E = (E B * + E MB )/2 as 
the mean of the energies of the states and 6 = Emb ~ E B * as their difference. Fixing an 
initial and final lattice timeslice ti and tf and summing over all possibilities for a single 
transition from B* to MB, which can occur at any intermediate time-slices between ti and 
tf, it follows that 


(B*, t f | MB, ti) = (B* | | MB) = (B* \ T nfi \ MB) 

Ufi-l 

— e -(E-8/2)t„ ^* | f | Mjg'j e -(E+5/2)(At fi -t n -a) _|_ 

^ +... (19) 


n =0 

= ax 


sinh(h At/j/2) Ea tfi 


sinh(a<5/2) 


with A tfi = tf — ti = anfi. In Eq. (19) the ellipsis denotes contributions of higher orders 
in the matrix elements (B* \ T | MB) and (MB \ T | B*), which are at least quadratic in the 
time separation A tfi. As a consequence, by extracting the term linear in Atfi we get the 
transfer matrix element 


1 MB) ~ —a(B* | H | MB). 


( 20 ) 


Given the time-dependence of the overlap in Eq. (20), we use two different fit ansatze 
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given by 


fi(t) = C 0 + Cl sinh (c 2 1/(2a)) /(c 2 /2), 

3 


t 


fl(t) — Co + Cl —b c 2 I — + 


( 21 ) 

( 22 ) 


In both cases, we are primarily interested in the parameter c±. Although in principle we 


could take into account the next terms denoted by the ellipsis in Eq. (22), in practice we 
will not need to go beyond c 2 to obtain a good fit to the available lattice data with their 
present accuracy. The parameter Co allows for an offset at t — 0, which can originate from 
lattice artifacts or contributions from excited states giving an overlap at zero time, i.e. with 


no insertion of the transfer matrix. Given that Eq. (19) contains a lattice version of the 


energy Dirac-h function for the finite temporal lattice extent, Eq. (21) allows for fitting the 
data taking into account a non-zero energy gap. Besides taking into account a finite energy 
gap, which will be the only contribution to it if the transition happens once in the path 
integral, the c 2 -term may also effectively be including next-to-leading order contributions 
arising from overlaps from other intermediate states. 

The value of the parameter c 2 extracted from fitting fi and the one extracted from fitting 
/ 2 can be related at order t 3 after expanding /j. We indeed find that these two values of 
the c 2 parameter show a strong correlation, as expected. Yet in some cases a different 
value is extracted and hence it is appropriate to use both ansatze to study the systematic 
uncertainties in the fitting of the lattice QCD data. 

As alluded to in Eq. (17), the interpretation of the overlap of lattice states and interpo¬ 
lating fields acting on the vacuum state involves the spinors ub* and ub as follows 

X] |(1 + 7o) (0 | Jg{t, x) e iP ' s | B*,P,s = 3/2, s 3 > = Z B *(P 2 ) VSpp, u^(P, s 3 ) 

X 

(23) 

Y j(l + T») <° Mt, V) J%{t, 2) 1 MB, k, s = 1/2, s 3 ) = Z„(P) Z B (k 2 ) V Sgg, u%(-k, s 3 ). 

V 

(24) 

We denote by s the spin quantum number of the baryon helds, sb * = 3/2 and sb = 1/2 
(having sm = 0 hxed for the pseudoscalar meson) and by s 3 its projection to a specific axis. 


The definition in Eq. (24) assumes that the total linear momentum of the meson-baryon 
state is zero. 
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While the ratio R^b i s constructed such that the numerical factors Z b */ m /b cancel, the 
spinors remain in the fraction and are combined to spin sums in the numerator and 

T, b * and T, b in the denominator via summation over the third spin component. To that 


end, we parametrize the slope c\ = ax in Eq. (21), (22) as follows 


Ci _ y- M(P = 0, |fc|,q~3,r 3 ) y ^mb(P = a 3 .r 3 ) 

ff3iT3 VN b * N mb ,/s B ,(P = 0)S B (|fc|) 


(25) 


'%b{P = o,|fc|^ 3 ) = J 


£ si § n 0') \ (! + 7o)a/J lifl. (-P = 0, cr 3 ) -U B 

° j=±l,±2,±3 




Cj; T3) ; 


where we use the notation e_j = — e} for j = 1, 2, 3 for brevity. The volume factor V stems 
from the lattice Kronecker —5 in momentum space for the total linear momentum P. We 
note that by our construction of the correlators, which are not summed over the source 
locations, but fulfill momentum conservation, we effectively have to insert this factor by 
hand for correct normalization of the ratio. N b * and Nmb denote the normalization of the 
states and we use the standard continuum-like normalization of on-shell states 


(B*,P, s 3 | £>*, P', S3) = Eb < P ' 2 ) v 5ss, 8 


771 B * 


PP> "S3S 3 


E B {kl) 


(MB, k M , k B , s 3 | MB, k' M , k' B , s' 3 ) — 2 E M (k 2 M ) V 8k M %' M x ~N~~ ^ ^s 3 s!, 

/ / i 


N b * — V , Nmb — 2 V 


2 E M (k 2 )E B (k 2 ) 


m B 


(26) 


In the last line of equation (26) we specialized to the case at hand with P = 0 and 
kM — ~k B = k. 


To sum the spinors in the numerator we likewise parametrize the matrix element according 
to leading order effective field theory, 

B* 

M(P B *, k M , k B ,a 3 ,r 3 ) = Cqg ~ AIB (Pb*, 03) kM^ u B (k B , r 3 ) . (27) 

ZTfl q 

Ccg = Cqg(Ib* j I B * I 0, 0; I B , I B ) is the Clebsch-Gordan coefficient for the coupling of the 
isospin of M and B to match that of B*. 


With Eq. (27) we can then use the standard spin sums for spin-1/2 fermions and the 
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B * In*, 4* 

M IMi I'm 

B I B , I 3 b 

Ccg| 

A++ 3/2, +3/2 

TT + 1 , +1 

N + 1/2, +1/2 

1 

£*+ 1, +1 

vr + 1 , +1 

A 0, 0 

1 

£*+ 1, +1 

7T+ 1, +1 

E° 1, 0 


1/2, -1/2 

TT~ 1 , —1 

H° 1/2, +1/2 



TABLE I. We give the isospin quantum numbers of the decuplet (first column) and the two-particle 
decay channel consisting of a meson with quantum number given in the second column and a spin- 
1/2 baryon with quantum numbers given in the third column. In the last column we give the 
absolute value of the isospin factor Ccg- 


Rarita-Schwinger field[36, 37] 


Eb(|*|) = 


E B (k 2 ) + m B 


2m 


B 


Y, b *(P — 0 ) — - 


Sf; B (P = 0, \k\) = S fl (|Al|) x E b .(P = 0) = ^ 


1 E B (k 2 ) + m B 


We thus write the coupling as 


B * _ a/ N b * Nmb 2 m B I 1 E B (k 2 ) + m B 
9mb ~ c i 


m B 


- 1/2 


(28) 


VC, 


CG 


\k\ 


m B 


(29) 


With this expression we can go back and rewrite the matrix element in terms of the extracted 
slope ci, 

,b* \ 2 


\M.(P — 0, |A:| )| = Cq G 


9 mb \ 2 p E B {k 2 ) + m B _ 2 2 N b *N mb 
2 m B ) 3 m B 1 V 2 


(30) 


We note, that the expression in Eq. (30) gives the squared matrix element for the transition 
between a certain isospin state of B* and a certain product of isospin states for M and B , 
such that 


C CG = </ 


B* 


I 3 

I B * 


'M: 


f 3 


1 B 


II). 


( 31 ) 


In Table |T] we give the isospin values for the decuplet resonances and their decay channel. 
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C. Density of states 


To apply Fermi’s Golden Rule we need to estimate the density of states at the transition 
energy. For a free pion (pseudoscalar M) and a free baryon B in the center-of-mass frame 
with k — |jfc| the total energy is 

Ef(k) = E M (k ) + E B (k) = \Jm M + k 2 + \J nn B + k 2 . 

Furthermore, we assume an isotropic density in the volume L 3 , with the unit cell in mo¬ 
mentum space being of size 27 t/L. Up to momentum k we thus count Q(k) / (2n / L) 3 = 
Airk 3 /3/(2n/L) 3 states. Varying k we then have a density of states given by 

L 3 


dfl L 3 L 3 2 dk 

p(Et) = dE,W = W k IF, 


8tt 2 


k 


E f 


(32) 


D. Decay width to leading order 


Having the overlap x from the lattice correlator functions and using the density of states 
we can connect, to leading order in the effective theory, the decay width in the continuum 
to that on the lattice by suitable normalization. To that end we observe, that to leading 
order in the continuum effective held theory, we have 


r = 


m B * 


2s B * 1 E 


B* 


y^ \m.{s b *, sb)i 2 I 


dE f 


m B 


$B* i s B 


k(E f y dk(E f ) 

E B 2Em 2iB dEf 


2tt 5(Ef — Ej) . 

(33) 


Evaluating the ^-f un ctional in the center-of-mass frame with Ei — m B * = E B * we obtain 

1 


T = 2 tt 


= 2vr 


2 s B * T 1 


y: \m (s B *, ss)i‘ 


S B* j S B 


U 3 1 

77 P(Ei) 


Nmb N b * V 


2 c\ 


2 s B * + 1 
= 27t (\M.\ 2 ) p(Ei). 


p{Ei) 


(34) 


We note that the expression of Eq. (34) contains the sum over all final states (in particular 
all spin configurations of the held B) and the average over all initial spin states of the spin- 


3/2 baryon B*. In Eq. (34), the width calculated is independent of the normalization of 


states chosen at intermediate stages, as one would expect. The coupling in Eq. (29), on 
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the other hand, carries an explicit dependence on the normalization of states shown by the 


appearance of the factor y/N B * N MB . 

We would like to note that in a realistic setup the lattice Hilbert space is of high dimension 
and the lattice transfer matrix correspondingly large. The restriction to a two-dimensional 
subspace may still be justified, if the first excited states in the B* and MB channels are 
sufficiently higher in energy. As usual this would lead to an exponential suppression of 


contributions from such states as assumed in the ellipsis in Eq. (19). Only in this case the 
overlap is proportional to the time separation, receiving contributions from single transitions 
from initial to hnal state anywhere along the time axis. 


III. NUMERICAL RESULTS 

We analyze two ensembles: one for a hybrid action with domain wall valence quarks 
on a Nf — 2 + 1 staggered sea [38] and m n = 350 MeV and one for a unitary action with 
Nf = 2 + 1 domain wall quarks [32] and = 180 MeV. Subsequently we will use the 
labels “hybrid” and “unitary” to distinguish results obtained using these two sets of gauge 
configurations. Results for the A resonance for the hybrid calculation have been reported 
in Refs. [28j [29] and thus we do not discuss them in detail here. 


A. Simulation details 

For the hybrid setup we use an ensemble of Nf = 2 + 1 staggered fermion configurations 
with the light quark mass corresponding to a pion mass of 350 MeV and the strange quark 
mass fixed to its physical value. This MILC ensemble is labeled as MILC-2864- , m010m050 
[39]. As valence quarks we consider domain wall fermions with the light bare quark mass ad¬ 
justed to reproduce the lightest pion mass obtained using Nf = 2 + 1 staggered quarks [40]. 
The valence strange-quark mass was set using the Nf = 3 ensemble by requiring the va¬ 
lence pseudoscalar mass to be equal to the mass of the Goldstone boson constructed using 
staggered quarks [3D 02!- For the unitary setup we use an ensemble of gauge configurations 
generated by the RBC-UKQCD collaborations with Nf = 2 + 1 domain-wall fermions and 
the Iwasaki gauge-action labeled as RBCJ)lp75JL32T64-m045m001 [32] . The simulation 
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action 

L 3 x T m.ps/ MeV a/ fm L 5 

-^conf -^src 

OAPe/AaPE ^Gaussian/AGaussian 

hybrid 

28 3 x 64 350 0.124 16 

210 4 indep. 

2.0/20 4.0/50 

unitary 

32 s x 64 180 0.143 32 

254 4 coh. 

2.5/25 0.5625/70 


TABLE II. Simulation parameters. The second, third, fourth and fifth columns give the lattice 
size, pion mass, lattice spacing and size of the fifth-dimension for the two ensembles considered in 
this work. The sixth and seventh columns give the number of configurations used for the analysis 
and the number of source positions per configuration. The last two columns give the parameters 
for the APE and the Gaussian smearing used in the construction of the interpolating fields. 

parameters for both cases are given in Table {[TT]} . For the hybrid ensemble we perform four 
independent measurements on 210 gauge configurations. The source locations for these mea¬ 
surements are separated by T/4 in time direction and the spatial coordinates are randomly 
chosen across the spatial volume. In the case of the unitary ensemble we use 4 independent 
propagators, which are inserted coherently into a single sequential source. Upon subsequent 
inversion of the Dirac operator the latter gives rise to a superposition of four sequential 
propagators at distance T/4 in time direction and thus four coherent sets of contractions. 

We use source- and sink-smearing on all interpolating fields to improve the overlap of 
our interpolating fields with the ground state. The forward and sequential propagators are 
smeared using Gaussian smearing with the APE smeared gauge links entering in the hopping 
matrix of the Gaussian smearing function. The smearing parameters for both lattices are 
given in Table {[TT]} . 

Inversions of the Dirac operator have been performed using the packages QUDA [43, 44] 
for the hybrid calculation and Qlua [45] using Moebius-accelerated domain wall fermions for 
the unitary action [46] . 

In Fig. (3] we show the energies of the states that are relevant for our calculation. The 
energies for zero and one unit of momentum \q\ = 2^\/L are shown. We use a notation 
analogous to reference [25j giving the one-particle interpolating fields a subscript labeling 
their momentum, i.e. hq denotes the pion-state with zero momentum, 7Ti with one unit 
of momentum etc. Likewise Hi N\ is the pion-nucleon state, where each interpolating field 
is constructed with one unit of momentum, while keeping zero total momentum. With 
the label rr B, where B = N, A, E, 5 we denote the sum of the individual energies of the 
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pion and the octet baryon B. The individual energies are determined from the two-point 
correlators of each particle. 


A —>nN £* —> n K/tt £ 



FIG. 3. Energies of the states entering the study of the decays for the m n = 180 MeV unitary 
ensemble. The black arrows mark the transitions we consider in this work (cf. Eq. Q). 

For the unitary ensemble we observe a near degeneracy of energy levels for the A and the 
7r N scattering state as well as for the E* and the n A scattering state. On the other hand, 
a significant energy gap exists between the E* and the 7r E scattering state, and the 5* and 
the 7T S scattering state. The situation is qualitatively different for the hybrid calculation, 
for which the relevant spectrum is shown in Fig. |4]. We observe a larger energy gap for all 
transitions under consideration, which is roughly the same in all cases. 

This qualitative difference arises from the larger values of m,. In our approximation 
E M b = E m + Eb, so the gap is 5 = E B * — E M b = E B * — E M — E B . Considering Fig. (4] 
for the case A — 1 tN one observes that at q = 0, En + E n is significantly greater than E& 
so 6 is negative and far from threshold. Taking into account that E^ + E^ increases with 
q, and so does |<5|, the gap gets even bigger as compared to |5| at zero q. The same thing 
happens for E* and for 5*. In contrast, in Fig. (3) for the ensemble at 180 MeV pion mass, 
we see that the A is unstable and there is a chance 5 will pass through zero at the relevant 
q for the transition. Indeed we see that 5 is small and slightly positive. Similarly, for E*, 5 
is small and slightly negative. 
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A — > 7T iV 


E* —> n A/7r E 



FIG. 4. Energies of the states entering the study of the decays for the m n = 350 MeV hybrid 
ensemble. The black arrows mark the transitions we consider in this work (cf. Eq. Q). 

B. A —y ttN 

We first discuss the case of the A resonance. As shown in Fig. [3], the lattice kinematics 
produce a scattering tiN state that is approximately degenerate with the A mass, thus 
satisfying one of the conditions for the validity of the method. For the hybrid action, used 
in our study, the energies have a sizeable gap as shown in Fig. |4]. In Fig. |5] we show the 
ratio R^ n for both the unitary and the hybrid action. 

For the unitary case we compare two different ways to combine the available lattice 
data. In the first one, all individual factors of the two-point correlators entering the ratio 
are averaged before the ratio is built using the maximal set of lattice symmetries for the 
individual correlators. This way of combining data will benefit from possible cancellation 
of additive lattice artifacts in individual correlation functions. We refer to it as maximal 
averaging. In the second approach, we build the ratio for each data set given by the tuple 
(momentum / direction, forward and backward propagators, source location) and in the final 
step combine the individual estimates for the ratio. Since the correlator data within one and 
the same tuple is maximally correlated, such an average would benefit error cancellations 
due to statistical correlation. 

We note, that due to the coherent source method used with the unitary action, we must 
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unitary, max. average * 
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0.01 


0.00 
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FIG. 5. Ratio R^ N for the unitary and hybrid calculation; the detail plot shows the evolution of 
the statistical uncertainty with t/a. The errors are computed using the T-method |5Tj. The red 
crosses are obtained using maximal averaging to construct the ratio for the unitary action, while 
the blue triangles exploit maximal correlation. The black triangles show the ratio obtained with 
the hybrid action using maximal correlation combination of data. Data points for different curves 
have been displaced horizontally. 


keep the source-sink time separation sufficiently smaller than the distance between the source 
insertions, i.e. tf «- U <C T/A. The ratio R// n shown in Fig. (5] exhibits a time dependence 
that is consistent with the expected linear behavior for both hybrid and unitary action, as 
well as for both types of averages. An overall comparison of the two approaches used for 
constructing the ratio with the unitary action does not reveal any significant difference in 
the mean value or the statistical uncertainty of the data points where they are both defined. 
However, on closer examination, the maximally averaged approach produces data for larger 
time slices. This is due to the fact that the A and nucleon correlators are more accurately 
determined having thus a lower probability of becoming non-positive in the sampling part 
of the error estimate. We shall thus use maximal averaging to combine data in what follows. 

In Fig. Im we show the results from fitting R^ N using the two fit ansatze given in 


Eqs. (21),(22), for the unitary action. We observe that the linear fit ansatz labeled as 


type 2, which uses a correlated fit with the function f 2 (t) and two free parameters c 0 , Ci, 
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already leads to fits with a value for x 2 /dof below one (bottom panel). The fit value for the 
slope determined by C\ does not show significant variation when scanning the fit ranges from 
[2, 5] to [6, 9]. As shown in Fig. [5], the statistical uncertainty of the fit parameters increases 
with increasing the lower fit range to larger values as expected from the dependence of the 
statistical uncertainty on the fitted data. Moreover, using as a fit function fi(t) we do not 
observe a significant change for c\. Neither do we observe any significant dependence of the 
central value of C\ on the number of parameters. In fact, C2 is statistically consistent with 
zero for both ansatze /1 and / 2 as shown in the center panel of Fig. [6). We point out that 
we include different labels on the left and right y -axes to show the values extracted using 
the fitting functions /1 and / 2 , respectively, because the order of magnitude of c 2 differs. 



fit range 

FIG. 6. We show the values of the parameters ci (upper panel), c 2 (central panel) and \' 2 /dof 
for the fit (lower panel) using three different ansatze to fit R^ N - The symbols represent from left 
to right: type 1 fits using fi{t) with 3 fit parameters (red squares), type 2 using / 2 (f) with 2 fit 
parameters (blue circles) and type 3 using / 2 with 3 fit parameters (black diamonds). The x-label 
gives the fit interval as a column t m i n atop f max in lattice units; in the case of c 2 we show results 
using type 1 with the labels to the left y-axis and type 2 with respect to right y-axis. 

As indicated in Fig. [6], we perform a large number of fits with different fit types using 
the fit ansatze / 12 (f) and different time ranges for the estimate of the slope and the energies, 
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for zero and one unit of momentum that enter the calculation. We note that with the A 
at rest, its mass does not enter the calculation of the coupling constant or the width. This 
result relies on energy conservation and generalizes to all decays studied here. 

Since many fits yield an acceptable value for y 2 /dof, we combine the different analyses 
with an appropriate weight to extract a mean value and meaningful estimates for the statis¬ 
tical and systematic uncertainty to account for the varying goodness of the fits and precision 
of the estimates from them. We consider the distribution of the results from each individual 
fit and associate a weight to it as follows: 

w(q, r ) = n w (^) 

A 

w(A ) = (1 — 2 |0.5 — p A \) 2 x var(A) -1 . (35) 


Here, pa denotes the p-value for the fit of quantity A, 


Pa I / x2 ,dofP0 dX , 


x 2 a 


where Xa * s the observed value of \ 2 for the fit of quantity A and f x 2 ;do f the density function 
for the x 2 distribution for dof degrees of freedom. A runs over all the quantities that have 
been derived from a fit and enter the calculation of g (or T), i.e. the slope parameter ci, the 


meson mass mM and baryon masses ms- The definition in Eq. (35) gives a higher weight to 


fits with a p-value close to 0.5, such that there is equal probability of Ending results above 
and below the observed fit value, and to those fits with smaller variance of the fit result. We 
then take a weighted average from the distribution as the mean value, 


9 = 


9iw(9i ) x w{9i)\ , 


where the sum runs over all fits labeled by index i. The statistical uncertainty is calculated 
from the variance of the bootstrap samples for the weighted mean, 

d C/stat = V Var (9) ■ 


Finally, the systematic uncertainty is estimated from the variance of the weighted distri¬ 
bution of the set {gi}- we form a histogram, where each g t gives a count proportional to w(gi) 
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to the corresponding bin. The square root of the variance derived from this distribution gives 
the systematic error 5 g sys . We then quote our results as 

9 = 9 9sta.t) (<5 9sys) ■ 


We proceed in the same way for the evaluation of the width T. 


we arrive at the values given in Eqs. (36) and (37). 


Following this procedure, 


9nN (unitary) = 23.7 (0.7) (1.1) 

(36) 

aV^ N (unitary) = 0.0868 (57) (32) 

(37) 


C. E* -> vrA 

For the decay of the E*, we follow the same approach as in the case of the A decay. In 
Fig. [7] we show the ratio for both the unitary and hybrid cases and Figs. |8) and (9] 
display the behavior of the three different fits when varying the fit ranges. For the unitary 
action, the energy levels for E* and n A are still close and we find acceptable linear fits 
already starting at t rn \ ri /a = 2. For the hybrid case we find, that starting with t m i n /a = 7 we 
find a time independent value for the slope even for the linear fit. The estimates for the slope 
from the cubic and hyperbolic sine fit are mutually consistent even before that. However, 
the x 2 /dof for all fit versions is acceptable starting as early as t miQ /a = 5. A qualitative 
difference between the unitary and the hybrid ensembles becomes apparent when examining 
the ratio . When attempting a linear fit to extract the slope, with the hybrid action, 
the central value for c\ rises systematically, when the lower end of the fit window is moved 
towards larger time-slices. This would be expected for a significant energy gap between the 
state excited by the E* and the 7r E interpolating fields. The upward curvature then shows 
that E n a > E^*. 


D. E*-^7tE 


For the transition E* —> 7tE we show the results for the ratio in Fig. 10 and the 


results for the parameters c 1 and C2 from a variety of our fits in Figs. 11 for the unitary 


and 12 for the hybrid calculation. 
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FIG. 7. Ratio R^+\ for the unitary and hybrid calculations. The notation is the same as that in 
Fig. [5] 
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FIG. 8. Fit range dependence of the fit parameters ci, C2 for three different fits of ratio for 

the unitary action. The notation is the same as that of Fig. [6j 
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FIG. 9. Same as Fig. [8j, but for the hybrid action. 


E. —> irE 


Finally we present the results for the transition S* —>■ 7rS in an analogous manner in 


Figs. 13 , 14 and 15 


We gather our results for the coupling and widths in Tables {|V|} and {VI} below. To 
allow for an easy comparison we convert the decay widths to physical units using the values 
for the lattice spacing given in Table ©• 


The results for the process A -H- 7riV with the hybrid action differ slightly from our 
previous investigation [28, [29], since we updated them using the weighted average for the 
distribution of fits. 


Utilizing the expressions of Eqs. (29) and (30) we estimate the coupling, which is in¬ 
dependent of the isospin combination of in and out state, while the width is for specific 
combinations of in and out states. For this reason, in the table we distinguish explicitly the 
isospin dependence of the width by giving the electromagnetic charges of the interpolating 
fields as superscripts. 
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FIG. 10. Ratio R ^+%o for the unitary and hybrid calculations.The notation is the same at that 
in Fig. [5j 
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FIG. 11. Fit range dependence of ci, C2 for three different fits of ratio R^+eo for the unitary 
action. The notation is the same as that in Fig. [6} 
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FIG. 12. Same as Fig. 


11 , but for the hybrid action. 


F. Improved precision with AMA for R^ N 


In order to assess the potential of using all-mode-averaging to improve the accuracy of our 
computations, we apply all-mode-averaging [3TLj on a subset of 89 (out of 254) configurations 
and specifically look at the case of A —y n N. In addition to the correlation functions, which 
had been obtained at high solver precision during the production of quark propagators, 
a corresponding data set at low solver precision was produced with 16 random shifts of 
the original spatial source position two-point correlation functions Ca-a, C n _ n and Cm-n 
for each of the four preset source time-slices independently. The measurements for Ca-ttAt 
are done coherently with a single inversion after inserting sequential sources at the four 
time-slices. 

We show a comparison of the estimates for the ratio R^ N and its statistical uncertainty 


in Figs. 16 and [17] . We find full consistency of the data for the ratio from both the AMA 
simulation and the original production run. Moreover, the uncertainty is reduced by a factor 
around two across the relevant time-slices 1 < t/a < 10. From an ideal scaling of the error 
we expect reduction of the statistical uncertainty by a factor of \J (89 x 64) / (254 x 4) « 2.4 
and the observed behavior is consistent with this expectation. 

We determine the coupling and width based on the AMA data set in the way previously 
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FIG. 13. Ratio for the unitary and hybrid calculations.The notation is the same as that 

in the Fig. [5} 



FIG. 14. Fit range dependence of c\ and C 2 from the fits of ratio B .~_^, 0 for the unitary calculation. 
The notation is the same as that in the Fig. [hj 
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FIG. 15. Same as Fig. 


14 , but for the hybrid calculation. 



r; A «F A 

ihr N aL ir N 

AMA 

24.2(0.3) (1.0) 0.0846(23) (22) 

HP, full 

23.7(0.7) (1.1) 0.0868(57) (45) 


TABLE III. The results on the A — irN coupling constant and width using AMA are shown in 
the first row, while results extracted without AMA are included in the second row. 


outlined and give the results for comparison in Table III 


The statistical errors show the 


expected improvement by a factor of approximately 2.4. 


IV. DISCUSSION 


In order to make a direct comparison with the experimental values, we provide in Table [TV| 
the data taken from the Particle Data Group [48] for the relevant baryon and meson masses, 
the full widths, branching ratios and relative momentum k of the asymptotic final meson 
and baryon state from the decay in the center-of-mass frame. The coupling constant is then 


derived according to the tree-level decay process using the expression in Eq. (33) and the 
experimental value of the width as an input. 
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FIG. 16. Comparison of results on the ratio R^ N with (red squares) and without using all-mode- 
averaging (blue circles and magenta triangles). 



FIG. 17. Comparison of the statistical uncertainties of with and without using all-mode- 
averaging; the detail plot shows the ratio of the uncertainty obtained by using AMA over the 
that from the original run with 254 configurations and 4 source time-slices per configuration; the 
horizontal line marks the ratio of uncertainties expected for ideal error scaling. 
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B * 

MB 

ms* / MeV 

uim / MeV m B / MeV 

r fu ii/MeV rfrVrfuu 

p/ MeV 

A 

nN 

1232(1) 

139.57018(35) 938.272013(23) 

118(2) 

1. 

227 

E* 

IT A 

1382.80 (35) 

1115.683(6) 

36.0(7) 

0.870(15) 

205 

E* 

7tE 


1192.642(24) 


0.117(15) 

120 

77* 

7r S 

1535.0(6) 

1314.86(20) 

9.9 (1.9) 

1. 

158 


TABLE IV. The physical values of masses mg*, mg and rriM, full widths Tf u n, branching ratios 
T fjB and relative momentum p for the 2-hadron state for the resonances studied in this work as 
given by the Particle Data Group [28]. 


process 

unitary hybrid 

PDG 

A ++ 7T+ N + 

23.7(0.7) (1.1) 26.7(0.6) (1.4) 

29.4(0.3) 

E* + O 7T+ A 

18.5(0.3) (0.5) 23.2(0.6) (0.8) 

20.4(0.3) 

S* + 7 r+ E° 

16.1 (0.3) (1.9) 19.0(0.7) (2.9) 

17.3(1.1) 

E*" W E° 

21.0(0.3) (0.3) 25.6(0.6) (4.3) 

19.4(1.9) 


TABLE V. Results for the couplings gfj B - F° r eac h decay process given in the first column, we give 
the coupling constant gf[ B for the unitary DWF ensemble with m n = 180 MeV (second column) 
and hybrid ensemble with m n = 350 MeV (third column). The fourth column shows the value of 
the coupling at leading order effective field theory using input from the PDG. The uncertainties in 


brackets are statistical and systematic as given in subsection III B 


We compare these values for the coupling constants to the results of our calculation 
in Table ||V|}. The analogous comparison for the decay widths in physical units is shown 
in Table ©■ We would like to stress that although we show the results for the hybrid 
and unitary calculation side-by-side in the tables, one should be careful in drawing strong 
conlcusions since the conditions for the applicability of the transfer matrix method are 
fulfilled to different degree in the two cases. In particular, the energy matching is very 
different in the two cases and for a direct comparison one would need to have kinematics 
where the energy gap is similar. 

We find that our lattice QCD values for the couplings are in good agreement with the 
PDG-derived values for all decays for both the unitary and the hybrid action with a tendency 
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process 

unitary 

hybrid 

PDG 

A ++ 7T+ N + 

119.4 (7.9) (4.5) 238.5 (12.2) (16.2) 

118(2) 

E* + -H> 7T + A 

54.5(2.1) (1.3) 

143.9(7.4) (6.1) 

31.3(8) 

E*+ o 7r+ S° 

17.6(0.8) (2.1) 

58.3(3.4) (6.8) 

4.2(5) 

E*“ 7 r“ S° 

35.1 (1.1) (0.4) 

126.0(5.6) (18.5) 

9.9 (1.9) 


TABLE VI. Results for the decay widths in MeV. The meaning of the columns is analogous to 
Table {|V|} with the decay process (first column), results from the unitary DWF ensemble with m^ 
= 180 MeV (second column), results from the hybrid ensemble with m n = 350 MeV (third column) 
and the PDG value (fourth column). 


of higher values for the latter case. This observed level of agreement is remarkable, given 
that with the unitary and hybrid action we simulate at pion mass 180 MeV and 350 MeV, 
respectively, and on coarse lattices. 


For the width itself, on the other hand, we only find agreement for T^ N with the unitary 
action. This may be expected since it is only for this case that the energies of the states 
B* and M B are degenerate and therefore this case is the closest to the threshold situation 


where the conditions of our approach are best fulfilled. Table |IV| shows in the right-most 
column the momentum in the center-of-mass frame for the fields M and B for the individual 
decays. On the lattice, this momentum is of course fixed to k = 2^t/L in lattice units or 
/e(unitary) rs 270 MeV for the unitary action and fc(hybrid) m 357 MeV for the hybrid one. 
Thus, in addition to matching the energies of the resonance and the decay channel, one has 
another constraint, namely, a fixed center-of-mass momentum in these transition processes, 
which deviates from the physical situation by a process-dependent amount. In general, we 
have that with the hybrid action, the lattice momentum is 1.5 to 3 times larger than its 
value in the continuum infinite volume limit. With the unitary action the violation is less 
severe, and the closest to the physical situation is the one corresponding to decay of the A. 


Assuming a finite volume we can check that the density of states derived from the lattice 
values of the masses and the momentum approaches closer to the value of the density of 
states derived with their continuum counterparts when going from n N to tiA/ttY to 7rS. 
This is to be expected, since the strange quark mass is tuned closer to its physical value 
than the light quark mass and A/S and 5 are have strangeness —1 and —2, respectively. 
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The dependence of the coupling and decay width on the meson and baryon masses, 
momentum and the parameters of the lattice simulation show a large disparity reflected in 
the different levels of agreement in Tables {fV{f and ©• Partly this is explained by the 
additional condition of having to match the center of mass momentum for extracting the 
width in the decay process. One would need to study the dependence of the momentum 
further in order to understand the different level of agreement between the case of the 
coupling and that of the width. 


V. CONCLUSIONS AND OUTLOOK 

The coupling constants g^ N , and g^ are evaluated using two ensembles of 

dynamical fermion gauge configurations with pion mass 350 MeV and 180 MeV. In both 
cases, domain wall valence quarks are used. The gauge configurations for the ensemble with 
the heavier mass were produced using Nf = 2 + 1 staggered sea quarks and thus our analysis 
is done with a hybrid action, while those with the lighter pion mass were produced using 
Nf = 2 + 1 domain wall sea quarks so the action is unitary. The kinematical conditions are 
best satisfied for the unitary action for all four decays, with the A-decay being closest to the 
physical situation. Comparing the values of the coupling constants obtained for these two 
ensembles, we find that they are about 10% smaller for the ensemble with 180 MeV pions as 
compared to their values for the ensemble with 350 MeV pions. Given that the pion mass is 
about half as compared to the hybrid ensemble, we conclude that the pion mass dependence 
is rather weak and thus the values obtained using 180 MeV pions should be close to their 
values at the physical point, which is indeed what we observe. In order to extract the width, 
one needs to make further assumptions, some of which are not well satisfied. For example, 
the energy in the center of mass frame on the lattice is different from the one in the infinite 
volume limit. The case for which these energies best match is for the A-decay where indeed 
we find an agreement with the experimental value. This demonstrates that the methodology 
works when the physical kinematical conditions are approximately satisfied on the lattice. 

To explore the applicability of AMA on reducing the statistical uncertainty of the ratio 
and consecutively of the extracted slope, we consider the all-mode-averaging technique for 
the case of the transition A +>• n N. Adding further correlation functions with randomly 
shifted source positions at low precision for a subset of the gauge ensemble, we increase 
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the available statistics by a factor of approximately 5.6. The ideally expected reduction 
of the statistical uncertainty is thus by a factor of 2.4. We observe an improvement of 
approximately a factor of 2 on our final, derived quantities, which is satisfactory. The solver 
precision for the low-precision inversions of the domain-wall Dirac operator is tuned to a 
compromise value that, on the one hand, yields sufficiently high statistical correlation for 
the two-point functions for both high and low precision inversions to ensure a good scaling 
of the statistical uncertainty with the number of low-precision inversions, and on the other 
hand, to keep the ratio of cost for a low- to a high-precision propagator as small as possible, 
which in our case turns out to be 1:5. 

Fully exploiting the potential of further reduction of the uncertainty of the slope bears 
the interesting prospect of becoming sensitive to contributions from excited states and next- 
to-leading order terms. This would be of particular importance in a more comprehensive, 
combined analysis of several decay channels and vital for an attempt to tackle the quark- 
connected diagrams, the calculation of which is beyond the scope of this work. Notwith¬ 
standing these future prospects, our current analysis shows, that for the time being, the 
major source of systematic uncertainty stems from the lattice kinematical setup rather than 
statistics. 

Given the good agreement of our lattice QCD results with the experimental values for 
the coupling constants and for the width when the kinematical constraints are satisfied, we 
plan to apply the method to study other baryon decays such as the decay of baryons in the 
negative parity channel and decays of baryon of higher spin. 

In the future we are also planning to address some of the deficiencies of the method con¬ 
nected to the kinematical conditions by considering moving frames. The decays considered 
here can be the test-bed for these extensions. 
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